THE JEANS MECHANISM AND 
BULK-VISCOSITY EFFECTS 



Nakia Carlevaro a ' b and Giovanni Montani 6 ' c ' di e 

"■Department of Physics, Polo Scientifico - Universita degli Studi di Firenze, 
INFN - Section of Florence, Via G. Sansone, 1 (50019), Sesto Fiorentino (FI), Italy 

b ICR A - International Center for Relativistic Astrophysics, 
c/o Dep. of Physics - "Sapienza" Universita di Roma 

° Department of Physics - "Sapienza" Universita di Roma, Piazza A. Mora, 5 (00185), Rome, Italy 

d ENEA - C.R. Frascati (Department F.P.N.), Via Enrico Fermi, 45 (00044), Frascati (Rome), Italy 

e ICRANet - C. C. Pescara, Piazzale delta Repubblica, 10 (65100), Pescara, Italy 

nakia. carlevaro@icra. it montani@icra. it 



Abstract: In this paper we study the gravitational instability in presence 
of viscosity. In particular, the standard Jeans Mechanism is analyzed taking 
into account bulk-viscosity effects in the first-order Newtonian dynamics. We 
show how the perturbation evolution is damped by dissipative processes and 
the top-down fragmentation mechanism is suppressed for strong viscous effects. 
The critical value of the Jeans Mass remains unchanged also in presence of bulk 
viscosity. 



1 Introduction 

The crucial dichotomy between the isotropy of region at red-shift z rs ~ I0 3 and the extreme 
irregularity of the recent Universe, z rs <C 1, is at the ground of the interest in the study 
of the gravitational instability for the structure formation. The fundamental result of the 
cosmological density-perturbation analysis is the so-called Jeans Mass, which is the threshold 
value for the fluctuation masses to condense, generating a real structure. If masses greater 
than the Jeans one are addressed, density perturbations begin to diverge as function of time 
giving rise to the gravitational collapse [H Ej ■ In this work, we consider dissipative effects 
into the fluid dynamics, in the linear Newtonian regime. The starting point is the fluid 
motion equations on which one can develop a first-order perturbative theory by adding small 
fluctuations to the unperturbed background. In particular, we introduce in the first-order 
analysis the so-called bulk viscosity, described as a function of the Universe energy density p 
via a power-law: Q = z p s where z,s = const. We show how the perturbations are damped 
by viscosity and the top-down fragmentation mechanism is suppressed. The critical value of 
the Jeans Mass remains unchanged also in presence of bulk viscosity. In the the standard 



Jeans Mechanism [U O [4] , the unperturbed background is assumed to be characterized by a 
static and uniform solution of the fluid parameters. 

2 Motion Equations of Viscous Fluids 

In order to describe the Newtonian evolution of a fluid, we here want to introduce the Eulerian 
equations governing the fluid parameters: the density p, the local 3- velocity v (of components 
v a ) and the pressure p, in presence of a gravitational potential <f> (indices a and (3 run from 1 
to 3). Since we are interested to treat isotropic and homogeneous perturbative cosmological 
models, we can safely neglect the so-called first viscosity {shear viscosity) in the unperturbed 
dynamics. In fact, in such models there is no displacement of matter layers with respect to 
each other and this kind of viscosity represents the energy dissipation due to this particular 
effect. Indeed, in presence of small inhomogeneities, such affect should be taken into account, 
in principle. However, in this work, we are aimed at studying the behavior of scalar density 
perturbations. In this respect, volume changes of a given mass scale are essentially involved 
and, therefore, we concentrate our attention to the effects of the so-called second viscosity 
or bulk viscosity, only. In fact, such kind of viscosity is responsible for the non-equilibrium 
dynamics of matter compression and rarefaction and we expect that such a phenomenon 
is more relevant than friction among different layers. According to literature developments 
S 13 [lOl [11] , we now assume the bulk- viscosity coefficient ( as a function of the 
energy density p(t) expressed via a power-law of the form 

C = z p s , z,s — const. , (1) 

here z is a constant parameter which defines the intensity of viscous effects |12) . 

Viscous fluids are governed, in Newtonian regime, by the following set of equations |13j : 

dtp + V • (pv) = , (2a) 
pd t v + p(v- V)v + Vp - CV(V-v) + pV0 = O, (2b) 

V 2 - AnGp = , (2c) 

which are the Continuity Eq. (energy conservation), the Navier-Stokes Eq. (momentum 
conservation) and the Poisson Eq. for the gravitational field </>, respectively. The pressure 
and density are linked by the Eq. of State (EoS), i.e., p = p{p)'. in this picture, the sound 
speed is defined by the relation v"l — g p/sp. 

3 Analysis of the dissipative Jeans Mechanism 

The mass agglomerates are due to the gravitational instability: if density perturbations are 
generated in a certain volume, the gravitational forces act contracting this volume, allowing 
a gravitational collapse. The only forces which contrast such gravitational contraction are 
the pressure ones. The Jeans Mechanism [T] analyzes what are the conditions for which 
density perturbations become unstable to the gravitational collapse. Such a model is based 
on a Newtonian approach and the effects of the expanding Universe are neglected. The 
fundamental hypothesis is the static and uniform solution for the zeroth-order dynamics 

v = , Po — const. , p = const. , <p — const. (3) 

Of course, this assumption contradicts the gravitational equation, but we follow the original 
Jeans analysis imposing the so-called "Jeans swindle" [U [3] . We underline that our study will 
focus on Universe stages when the mean density is very small: in particular the recombination 
era, after decoupling. This way, the effects of bulk viscosity on the unperturbed dynamics 
can be consistently neglected in view of its phenomenological behavior (fT]). 
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3.1 Review of the Jeans Mechanism 



In the original Jeans Model, a perfect fluid background is assumed. After setting C = in 
eqs. (0), let us add small fluctuations to the unperturbed solution: p = p„ + Sp, p = p„ + Sp, 
<f> = (fi a + 5<f), and v = v + 5v. Substituting such expressions, one differential equation for 
the density perturbations can be derived: 

d?5p - v 2 V 2 Sp = AttG p Sp . (4) 

To study the properties of Sp, we now consider a plane-wave solutions of the form 

Sp(r,t) = ^e wt - tkT , (5) 

where u> and k (k = |k|) are the angular frequency and the wave number, respectively. This 
way, one can obtain the following dispersion relation 



uj 2 



= v 2 k 2 -4irGp . (6) 



In this scheme, two different regimes are present: if ui 2 > a pure time oscillatory behavior 
for density perturbations is obtained. While if to 2 < 0, the fluctuations exponentially grow 
in time, in the t — > oo asymptotic limit (i.e., we choose the negative imaginary part of the 
angular-frequency solution) and the gravitational collapse is addressed since also the density 
contrast S — s P/p diverges. The condition uj 2 = defines the so-called Jeans Scale Kj and 
the Jeans Mass Mj (which is the total mass in a sphere of radius R = ir/Kj). Such threshold 
quantities read 



3.2 Viscous effects 

Let us now analyze how bulk viscosity can affect the gravitational collapse dynamics and 
we recall that we are able to neglect such kind of viscosity in the unperturbed dynamics, 
which results to be described by the static and uniform solution ^ . We start by adding the 
usual small fluctuations to such a solution, i.e., Sp, Sp, 8(f), S~v. In treating bulk-viscosity 
perturbations, we use a standard Taylor expansion £ = Co + where £ = C(/°o)- Substituting 
all fluctuations in the system ([2]), one can obtain an unique equation for density perturbations, 
describing the dynamics of the gravitational collapse: 

Po d 2 5p ~ Po v 2 s V 2 Sp - Co V 2 d t 5p = AnGpl Sp . (8) 

Plane-waves solutions §E§ can be addressed, obtaining a generalized dispersion relation 

Po uj 2 ~ i Co k 2 uj + poiinGpo ~ v 2 k 2 ) = . (9) 

As in the standard Jean Model, the nature of the angular frequency is responsible of two 
different regimes. The dispersion relation has the solution 

oj = i(C a k 2 )/(2p )±V^ , u = -(k 4 C 2 )/(Ap 2 )+v 2 k 2 ~A7rGp , (10) 

thus we obtain the time exponential regime for to $C and a damped oscillatory regime for 
Q > 0. The equation uj — admits the solutions K\ and K2 which read 

*i,2 = (V2p v s )/(C„) (1 T y/l ~ (Kj{ )/(p v s ) 2 ) 1/2 , K 1 ,K 2 >0, K x < K 2 . (11) 
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The existence of such solutions gives rise to a first constraint for the viscosity coefficient: 
Co *S Cc = PoVs/Kj. An estimation in the recombination era0, after decoupling, yields to 
the value Cc = 7.38 • 10 4 gem -1 s _1 and confronting this threshold with usual viscosity (e.g., 
^Hydr. _ g ^ 10~ 7 g cmT 1 s _1 ), we can conclude that the range Co ^ Cc is the only of physical 
interest. Finally we obtain: Q ^ for k ^ K\, Ky. ^ k and Q > for Ki < k < K 2 . 
Let us now analyze the density-perturbation exponential solutions for u> ^ 0: 

5p~e wt , w = -(C„fc 2 )/(2 Pn ) T V Z ^. (12) 

To obtain the structure formation, the amplitude of such stationary waves must grow for 
increasing time. The exponential collapse for t — ► oo is addressed, choosing the (+) sign 
solution, only if w > 0, i.e., k < Kj with Kj < K\ < K 2 . As a result, we show how the 
structure formation occurs only if M > Mj, as in the standard Jeans Model. The viscous 
effects do not alter the threshold value of the Jeans Mass, but they change the perturbation 
evolution and the pure oscillatory behavior is lost in presence of dissipative effects. In par- 
ticular, we get two distinct decreasing regimes: for K\ < k < K 2 (i.e., u> > 0), we obtain a 
damped oscillatory evolution of perturbations: 5p ~ exp[— t(( k 2 ) / (2p )] cos (y/(5 t) , while, 
for Kj < k < K 1 and K 2 < k, density perturbations exponentially decrease as Sp ~ e wt , 
with w < 0, in the limit t — * 00. 

4 Implication for the top-down mechanism 

As shown above, since the pure oscillatory regime does not occurs, we deal with a decreasing 
exponential or a damped oscillatory evolution of perturbations. This allows to perform a 
qualitative analysis of the top-down fragmentation scheme [14] . i.e., the comparison between 
the evolution of two structures: one collapsing agglomerate with Mq 3> Mj and an internal 
non-collapsing sub-structure with Ms < Mj. If this picture is addressed, the sub-structure 
mass must be compared with a decreasing Jeans Mass since the latter is inversely proportional 
to the collapsing-agglomerate background mass. This way, as soon as such a Jeans Mass 
reaches the sub-structure one, the latter begins to condense implying the fragmentation. 
In the standard Jeans Model, this mechanism is always allowed since the amplitude for 
perturbations characterized by Ms < Mj remains constant in time. On the other hand, the 
presence of decreasing fluctuations in the viscous model, requires a discussion on the effective 
damping and an the efficacy of the top-down mechanism. Of course, such an analysis contrasts 
the hypothesis of a constant background density, but it can be useful to estimate the strength 
of the dissipative effects. We now study two cases for different values of the bulk-viscosity 
coefficient: Co <C 1 and Co > 1- I n this analysis, a perturbative validity-limit has to be set: we 
suppose s P/p ~ 0.01 as the limit of the model and we use the recombination era parameters 
(see footnote 1 ), in particular the initial time of the collapse is define as the beginning of the 
matter-dominated Universe, i.e., t = tuD = 1-39 • 10 13 s. 

In correspondence of a very small viscosity coefficient (Fig. I), we consider the Jeans Mass 
Mj = 10 6 M Q and a sub-structure wave-number Ks in the region K\ < Ks < K 2 , i.e., 
density perturbations evolve like damped oscillations. Fluctuations have to be imposed small 
at the initial time t , this way, we consider density contrasts (5q for the galaxy and 5s for 
the sub-structure) of 0(1CP 3 ). In this scheme, the galaxy starts to collapse and the validity 
limit is reached at t* — 6.25 • I0 13 . As a result, in Fig.l we can show how the sub-structure 
survives during the background collapse until the threshold time value t*. Thus, we can 
conclude that, if the viscous damping is sufficiently small, the galaxy formation occurs. 

"^The parameters are set as follows: the usual barotropic relation p — c 2 / p 1 ^ 1 is assumed and the constant p can 
be derived from the expression expression Mj ((VJ. Universe is dominated by matter and we can impose the values: 
Mj ~ 10 6 A/ o , 7 = 5/3, p c = 1.879 h 2 ■ 10" 29 g cm~ 3 , h = 0.7, z = 10 3 and p = p c z 3 = 0.92 • 10 -20 g cm~ 3 . 
Using these quantities one finds p — 9.034 ■ 10 — 7 gcm~ 3 , v s — 8.39 ■ 10° cms -1 and the threshold value £ c . 
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Figure 1: Case Co = 14 g cm 1 s 1 . Galaxy density 
contrast: <5 G - M G — 10 12 Mq - (dashed 
line). Sub-structure density contrast <5g - 
Ms — Mq - (normal line). 



Figure 2: Case Co = 10~ 5 gem- 1 s _1 . Galaxy den- 
sity contrast: <5 G - M G = 10 12 M - 
(dashed line). Sub-structure density con- 
trast 5s - Ms — 10 Mq - (normal line). 



Let us now discuss the case ( > 1 (Fig. 2) by changing the sub-structure mass, which is now 
Ms = Mq. Here, the damping effects is stronger. In fact, when the galaxy density contrast 
reaches the threshold value Sg = 0.1, we obtain 6s = 10 -5 . The top-down mechanism 
for structure formation results to be unfavored by the presence of strong viscous effects: 
the damping becomes very strong and the sub-structure vanishes during the agglomerate 
evolution. 

5 Concluding remarks 

The effects induced by the presence of bulk viscosity on the Jeans Mechanism have been 
analyzed in a perturbativc scheme. Viscosity has been introduced in the first-order dynamics 
via a power-low function of the energy density and the effects produced on the unperturbed 
dynamics have been consistently neglected in view of the phenomenological nature of such 
kind of viscosity. 

The main result has been to show how bulk viscosity damps the density contrast evolution 
maintaining unchanged the threshold value of the Jeans Mass. Such an effect suppresses 
the sub-structure formation in the top-down fragmentation mechanism. In particular, a new 
decreasing regime for perturbations has been found. The presence of such a behavior allowed 
the study of the top-down scheme for small and strong viscous effects. In the first case, the 
density-perturbation amplitude of a sub-structure remains substantially constant during the 
main structure collapse. On the other hand, if viscous effects are sufficiently strong, the 
sub-structure vanish in the linear perturbative regime, unfavoring the fragmentation. 
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